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MACDONALD-KOORNWINDER MOMENTS AND THE 
TWO-SPECIES EXCLUSION PROCESS 


SYLVIE CORTEEL AND LAUREN K. WILLIAMS 

Abstract. Introduced in the late 1960’s [MGP681 [Spi70| , the asymmetric exclu¬ 
sion process (ASEP) is an important model from statistical mechanics which de¬ 
scribes a system of interacting particles hopping left and right on a one-dinrensional 
lattice with open boundaries. It has been known for awhile that there is a tight con¬ 
nection between the partition function of the ASEP and moments of Askey-Wilson 
polynomials [USW041 ICW1U ICSSW121 , a family of orthogonal polynomials which 
are at the top of the hierarchy of classical orthogonal polynomials in one variable. 
On the other hand, Askey-Wilson polynomials can be viewed as a specialization of 
the multivariate Macdonald-Koornwinder polynomials (also known as Koornwinder 
polynomials), which in turn give rise to the Macdonald polynomials associated to 
any classical root system via a limit or specialization |vD95j . In light of the fact that 
Koornwinder polynomials generalize the Askey-Wilson polynomials, it is natural 
to ask whether one can find a particle model whose partition function is related to 
Koornwinder polynomials. In this article we answer this question affirmatively, by 
showing that Koornwinder moments at q = t are closely connected to the partition 
function for the two-species exclusion process. 
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1. Introduction 


Introduced in the late 1960’s [MGP681 !Spi70| , the asymmetric exclusion process 
(ASEP) is a model of interacting particles hopping left and right on a one-dimensional 
lattice of N sites. In the most general form of the ASEP with open boundaries, parti¬ 
cles may enter and exit at the left with probabilities a and 7 , and they may exit and 
enter at the right with probabilities /3 and 6. In the bulk, the probability of hopping 
left is q times the probability of hopping right. The ASEP is important in statistical 
mechanics because it is one of the simplest models which exhibits boundary-induced 
phase transitions. Moreover, the ASEP has been cited as a model for traffic flow and 
protein synthesis. 

It has been known since work of Uchiyama-Sasamoto-Wadati [ IJSW04] that there 
is a close connection between the partition function Zn of the ASEP, and Askey- 
Wilson polynomials, a family of orthogonal polynomials which are at the top of the 
hierarchy of classical orthogonal polynomials in one variable. Using their work, we 
showed in [CSSW12] that each Askey-Wilson moment equals a specialization of the 
fugacity partition function Z N (f) of the ASEP, see Theorem O In jCWTTi fCWl3] . 
we introduced some new combinatorial objects called staircase tableaux, and used 
them to completely describe the stationary distribution of the ASEP; in particular, 
the partition function -Zjv(£) can be written as a sum over staircase tableaux of 
size N. It follows that (up to a scalar factor), Z^{f) is a polynomial with positive 
coefficients, and that Askey-Wilson moments can be expressed in terms of staircase 
tableaux. 

Askey-Wilson polynomials can be viewed as a specialization of the multivariate 
Macdonald-Koornwinder polynomials, which are also known as Koornwinder poly¬ 
nomials [Koo92j . or Macdonald polynomials for the type BC root system [ Mac95j . 
(For brevity, we will henceforth call them Koornwinder polynomials.) These polyno¬ 
mials are particularly important because the Macdonald polynomials associated to 
any classical root system can be expressed as limits or special cases of Koornwinder 
polynomials lvD95l . Since Koornwinder polynomials generalize Askey-Wilson poly¬ 
nomials, and Askey-Wilson moments are closely connected to the ASEP, it is natural 
to ask if there is some particle model generalizing the ASEP whose partition func¬ 
tion is related to Koornwinder polynomials. This question was posed to us by Mark 
Haiman in 2007 [Halj . 

When q = t, Koornwinder polynomials can be expressed in terms of Askey-Wilson 
polynomials by means of a Schur-like determinantal formula. This fact, together with 
the connection between Askey-Wilson moments and the partition function Zn(0 
of the ASEP, led Eric Rains IRaij to suggest that we consider, for any partition 
A = (Ai, A 2 , • • •, A m ), the quantity 



( 1 ) 


The K x (ty S can be considered to be Koornwinder moments, see Sections 121 and PT7T1 
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In this paper we will answer Haiman’s question affirmatively by demonstrating 
that there is a close connection between Koornwinder moments and the partition 
function of the two-species ASEP. The two-species ASEP is a generalization of the 
ASEP which involves two different kinds of particles, “heavy” and “light”. Both 
types of particles can hop left and right in the lattice (heavy and light particles 
interact exactly as do particles and holes in the usual ASEP), but only the heavy 
particles can enter and exit the lattice at the left and right boundary. So in particular, 
the number of light particles is conserved. When there are no light particles, the 
two-species ASEP reduces to the usual ASEP. The main result of this paper is an 
interpretation of the “complete homogeneous” Koornwinder moments A'( m ,o,o,...,o)(0 
in terms of the partition function of the two-species ASEP. That is, we prove the 
following. 


Theorem 1 . 1 . The Koornwinder moment K(jv-r,o,o,...,o)(0 (where there are precisely 
r 0’s in the partition) is proportional to the fugacity partition function Z N ^ r {f) for the 
two-species ASEP on a lattice of N sites with r “light” particles. More specifically, 

K(N-r, o,o,...,o)(£) = ^ _ qy^N, r (£)- 


(This formula makes sense when q ^ 1. For the q 


1 version, see Lemma \8J\ ) 


We also prove a Jacobi-Trudi-type formula which expresses any arbitrary Koorn¬ 
winder moment K\ as a determinant in the complete homogeneous Koornwinder 
moments, see Corollary 15.21 

Since the partition function is a polynomial with positive coefficients |CWlll 
ICW13] . Rains conjectured that the Koornwinder moments K\ are also polynomi¬ 
als with positive coefficients (up to a simple scalar factor) jRaij . Note that by the 
probabilistic interpretation of Theorem 11.11 it follows that AT(jv_r,o,o,...,o) is positive 
whenever we specialize the parameters a, (5,^,5, and q to be positive numbers be¬ 
tween 0 and 1. Moreover, in a sequel to this paper [CMW ]. we will give a tableaux 
formula for Z^.r, which implies that both Z^, r and K(N-r,o,o,...,o) are (up to a scalar 
factor) polynomials in the parameters with positive coefficients. 

It is worth noting that there has been some recent work which may be related 
to ours. In (BC14j . Borodin and Corwin introduced what they call Macdonald pro¬ 
cesses, which are probability measures on sequences of partitions defined in terms 
of nonnegative specializations of the (type A) Macdonald symmetric functions, and 
showed that they are related to various interacting particle systems. None of the par¬ 
ticle systems they discuss is the ASEP (or two-species ASEP) with open boundaries, 
but perhaps there exists some more general Koornwinder processes which would be 
connected to the ASEP. 

Additionally, while we were writing up this paper, we learned about some work 
of Cantini [Canl5j , which contains a result similar to our Theorem 11.11 However, 
Cantini uses the partition (l iV_r ,0 r ), and his techniques are completely different 
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from ours; e.g. he uses the affine Hecke algebra of type Cn as opposed to the Matrix 
Ansatz and the combinatorics of Motzkin paths. 

The structure of this paper is as follows. In Section [2j we give an introduction 
to Askey-Wilson polynomials, and Koornwinder polynomials, and Koornwinder mo¬ 
ments at q = t. In Section [3] we define the asymmetric exclusion process and the 
two-species exclusion process. We also explain the Matrix Ansatz, which is a pow¬ 
erful tool for analyzing the stationary distribution of these models. In Section U 
we give an interpretation of Koornwinder moments in terms of weighted Motzkin 
paths. In Section [5] we state and prove a Jacobi-Trudi type formula for Koornwinder 
moments. In Section [6] we prove a connection between partial Motzkin paths and the 
two-species exclusion process. In Section [3 we complete the proof of our main result. 
Finally in Section [8] we show that when we specialize £ = q — 1, the Koornwinder 
moments K\ have a beautiful multiplicative formula in terms of the hook lengths 
of the corresponding partition. This formula provides evidence for the positivity 
conjecture regarding Koornwinder moments. 

Acknowledgments: We would like to thank Mark Haiman, who pointed out to 
us in 2007 that Koornwinder polynomials generalize the Askey-Wilson polynomials, 
and asked us if we could make a connection between Koornwinder polynomials and 
some generalization of the ASEP. We would also like to thank Eric Rains, who 
suggested that we look at Koornwinder polynomials for q = t, and in particular at 
the quantities defined by (jT]) . Finally we would like to thank Jennifer Morse and 
Dennis Stanton for useful comments. 


2. Askey-Wilson polynomials and Koornwinder polynomials at q = t 


The Askey-Wilson polynomials are orthogonal polynomials with five free parame¬ 
ters (a, b, c, d, q). They reside at the top of the hierarchy of the one-variable orthogo¬ 
nal polynomial family in the Askey scheme [ AW85 . GR041IKLSIO] . In this section we 
define Askey-Wilson polynomials, following the exposition of AW 851 and [ USW04] . 
as well as Askey-Wilson moments. We will then define Koornwinder polynomials 
and their moments for q = t. 

The g-shifted factorial is defined by 

s n —1 

(ai, a 2 , • • • , a s ; q) n = JJ JJ (1 - a r q k ), 

r =1 k =0 

and the basic hypergeometric function is given by 


CXj ^ Qjj* 

h, ■ ■ ■ ,b s ' 


q,z 


E ( a l, ' ' ' , a ri q)k // -|\fc fc(fc-l)/2\ 

fa,---, J J 


k =0 


1 +s—rk 

/C • 


The Askey-Wilson polynomial P n (x) 
P n (x) = a~ n (ab, ac, ad; q) n 4 0 3 


= P n (x: a, b, c, d\q ) is explicitly defined by 
q ~ n , q n ~ 1 abcd, ae l6 , ae~ l6 

, , \q,q , 

ab, ac, ad 
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with x = cos 9 for n G Z + := {0,1, 2, • • • }. It satisfies the three-term recurrence 
A n P n+1 (x ) + B n P n (x) + C n P n _i(x) = 2 xP n (x), 
with P 0 (x) = 1 and P-i(x) = 0, where 
1 — q n ~ 1 abcd 

(1 — q 2n ~ 1 abcd)( 1 — q 2n abcd ) ’ 

n-i 

= (i - ^,mi - ,-Ai) 1(1++ “ Ws,) - 9 "" I(1+ '' )aW(s+9S ' )] ' 
(1 — g”)(l — g n_ 1 a&)(l — g" _ 1 ac)(l — g” _ 1 a<i)(l — g n_ 1 f>c)(l — q n ~ l bd)( 1 — q n ~ l cd ) 

(1 — g 2 r,_ 2 a&c<i)(l — g 2 n_ 1 a&cof) 


and s = a + 6 + c + d, s' = a 1 + b 1 + c 1 + d 1 


Remark 2.1. It is obvious from the three-term recurrence that the polynomials P n (x) 
are symmetric in a, b, c and d. 


For |a|, | 6 |, |c|, |d| < 1, using z = e , the orthogonality is expressed by 


c 


dz 

Pniz 


-w 


z + z 


-1 


z + z 


-1 


Pr, 


z + z 


-1 


h n ~ 

On 

hc\ 


where the integral contour C is a closed path which encloses the poles at z — aq k , bq k , 
cq k , dq k (k G Z + ) and excludes the poles at z — (ag fc ) _1 , (&g fc )” 1 , (cg fc ) _1 , (dg fc ) _1 
(k G Z + ), and where 


( = _ (e 2i6> , e~ 2ie ] g)oo _ 

(ae l9 , ae~ l9 , be l9 , be ~ %9 , ce l9 , ce~ l9 : de l9 , de~ %9 ] q , ) 0O ’ 
h n (1 — q n ~ 1 abcd ) (g, ab, ac, ad, be, bd, cd ; g) n 
ho (1 — q 2 ™- 1 abed) (abed] q) n 

, __ (abed; q ) 00 _ 

0 (g, ab, ac, ad, be, bd, cd] q)^ 

(In the other parameter region, the orthogonality is continued analytically.) 


Remark 2.2. Note that our definition of the weight function above differs slightly 
from the definition given in [AW85 ] ; the weight function in [AW85] did not have the 
ho in the denominator. Our convention simplifies some of the formulas to come. 


Definition 2.3. The moments of the (weight function of the) Askey-Wilson poly¬ 
nomials - which we sometimes refer to as simply the Askey-Wilson moments - are 
defined by 


hk = fJ>k(a,b, c, d\q) = 


dz 
47 t iz 


-w 


z + z 


-1 


z + z 


-1 
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Definition 2.4. Let z = (zi, ..., z m ), X = (Ai,..., A m ) be a partition, and a, b, c, d, q, t 
be generic complex parameters. The Koornwinder polynomials P\(z] a,b, c, d\q,t) 
are multivariate orthogonal polynomials which are the type BC-case of Macdonald 
polynomials. More specifically, P\(z] a,b,c, d\q,t) is the unique Laurent polynomial 
invariant under permutation and inversion of variables, with leading monomial z x , 
and orthogonal with respect to the Koornwinder density 


n 


l<z<j<m 


(zjZj, Zj/Zj, Zj/Zj, 1/zjZf, g) o, 
(tZiZjfiZi/zj, tZj/zi, t/ziZj■ q ) 


n 


{ZiA/zhv) c 


°o 1<i<m ( aZi , a/zi, bzi, b/z h cz u c/z h dz h d/zp, q) c 


on the unit torus \zi\ = \z 2 \ = ■ ■ ■ = \z m \ = 1, where the parameters satisfy 

Ml, Ml, Ml, Ml, Ml, Ml < I- 

At q = t, we have 


P A (z; a, b, c, d\q , q) = const • 


detp m _ j+Aj (^; a, b, c, d\g))™ j=1 
det (pm-fizi, a , b, c, d\q))™ j=l 


where the pi’s are the Askey- Wilson polynomials. 

Note that when q = t, the Koornwinder density becomes 


( 2 ) 


n 

l<i<j<m 


Zi/Zj) (1 - Zj/Zi)(l - 1/ZiZj) JJ w 

1 <i<m 



) 


where w denotes the Askey-Wilson density. 


Remark 2.5. When q = t, Koornwinder polynomials are sometimes called Mac¬ 
donald’s 9th variation of Schur functions associated with Askey-Wilson polynomials 
[NNSY01] . 

For Askey-Wilson polynomials, the kth moment pk is defined to be the integral of 
x k (here x = —) with respect to the Askey-Wilson density. For the multivariate 

Koornwinder polynomials, there are several ways that we conld define moments. One 
way would be to integrate a monomial in x±,... ,x m (here we set ay = " a+ ~ i ) with 
respect to the Koornwinder density (J2]). Following a suggestion of Eric Rains }Rai| , 
we will instead define our Koornwinder moments by integrating Schur functions in 
the ay’s with respect to (l2l) . 

Definition 2.6. Let /^(/(xi,..., x m )\ a, b, c, d; q, q) denote the result of integrating 
the function /(ay,... ,x m ) with respect to the Koornwinder density (J2]). We define 
the Koornwinder moment 

M a = I k (s x (x i,..., x m ] a, b , c, d; q, q)). 

As Rains pointed out to us, these Koornwinder moments have a determinantal 
formula [Ra il. 

Lemma 2.7. We have that 

M _ det (p,\ i+m -i +m -j ) ™j = i 

X det(/X 2 m-i-i)ij=l 
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where /ik is an Askey-Wilson moment. 


Proof. Note that when q — t, the Koornwinder density can be written (up to a scalar 
factor of as 

(3) Yl ( x i~ x i) 2 II 

l<i<j<m 1 <i<m 

where x t = ~ ,+ ^ . 

Recall that the classical dehnition of the Schur polynomials says that 


S\ (^i , , Xm) 


det(^ + — )% =1 

det(x; n ~% =1 


so that 


TT (Xi-Xj) 2 = det (x Xi+m *) det (x™ ‘). 

l<i<j<m 

Therefore we have that 

M\ oc [ det (x x f +m ~ l ) det (a;” 1- *) TT wixf) 

J 1 <i,j<m J 1 <i,j<m J ' i-i- 

Ki<m 


oc det 


xK+m-ixm-jw^x) 


det /^Ai+m—z+m —j • 


ffere we obtained the second line by applying the integral version of the Cauchy- 
Binet formula }dB55i IAnd83j . and we obtained the third line by using the definition 
of Askey-Wilson moments. 

Since the constant of proportionality is independent of A, we can recover it by 
setting A = 0, and using / fc (l) = 1. This gives us 


I K {s x {xi, ...,x m )\a,b, c, d\ q, q) 


det (/i 


\i+m—i+m—j 




m 

i,j =1 


det (/i 2 m-i-j )ij= 1 


□ 


3. The asymmetric exclusion process and the two-species exclusion 

PROCESS 

We start by defining the asymmetric exclusion process (ASEP) with open bound¬ 
aries. We will then define the two-species exclusion process, which generalizes the 
usual ASEP. Finally we will explain the Matrix Ansatz, which has been an important 
tool for analyzing the stationary distribution of the ASEP. 
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3.1. The asymmetric exclusion process (ASEP). 

Definition 3.1. Let a, j3, 7 , 5, q, andu be constants such that 0 <a<l, 0</3<l, 
0 < 7 < 1, 0 < 5 < 1, 0 < g < 1, and 0 < u < 1. Let Bn be the set of all 2 N 
words in the language {o, •}*. The ASEP is the Markov chain on B N with transition 
probabilities: 

• If X = A • o B and Y — A o »B then Px,y — 777 (particle hops right) and 
Py,x = 777 (particle hops left). 

• If X = o B and Y = »B then Px,y = 777 (particle enters from the left). 

• If X = B* and Y = Bo then Px,y = 771 (P ar ^c/e ex ^ s t° the right). 

• If X = *B and Y = oB then Px,y = 777 (particle exits to the left). 

• If X = Bo and Y = B» then Px,y — 777 (particle enters from the right). 

• Otherwise Px,y = 0 for Y ^ X and Px,x = 1 — x^y ^x,y- 

See Figure |T] for an illustration of the four states, with transition probabilities, for 
the case N = 2. The probabilities on the loops are determined by the fact that the 
sum of the probabilities on all outgoing arrows from a given state must be 1 . 


o 



13 


FIGURE 1. The state diagram of the ASEP for N = 2 


I 11 the long time limit, the system reaches a steady state where all the proba¬ 
bilities P n (ri, t 2 , ..., Tn) of finding the system in configurations ( 7 , r 2 ,..., Tn) are 
stationary. More specifically, the stationary distribution is the unique (up to scaling) 
eigenvector of the transition matrix of the Markov chain with eigenvalue 1. 

3.2. The two-species exclusion process. The two-species ASEP is a generaliza¬ 
tion of the ASEP which involves two kinds of particles, heavy and light. We will 
denote a heavy particle by a 2 and a light particle by a 1. We will also denote a 
hole (or the absence of a particle) by a 0. In the two-species ASEP, heavy particles 
behave exactly as do particles in the usual ASEP: they can hop in and out of the 
boundary, and they can hop left and right in the lattice (swapping places with a hole 
or with a light particle). Light particles cannot hop in and out of the boundary, but 
they may hop left and right in the lattice (swapping places with a hole or with a 
light particle). Therefore the number r of light particles is conserved. 

Definition 3.2. Let a, j3, 7 , 5, q, andu be constants such that 0 <a<l,0</3<l, 
0 < 7 < 1, 0 < S < 1, 0 < q < 1, and 0 < u < 1. Let Bx, r be the set of all words in 
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{0,1, 2}^ which contain precisely r 1 ’s; note that \B N)7 .\ = )2 JV r . The two-species 
ASEP is the Markov chain on Bx, r with transition probabilities: 

• If X = A21B and Y = A12B, or if X = A20B and Y = A02B, or if 

X = A10B and Y = A01B, then Px,y = 7777 and Py x = 7777- 

• If X = OB and Y = 2 B then Px,y = 77 ^ 77 - 

• If X = B2 and Y = BO then Px,y = 7777 - 

• If X = 2B and Y = OB then Px,y = 7777 ' 

• If X = BO and Y = B2 then Px,y = 7777 ■ 

• Otherwise Px,y = 0 for Y ^ X and Px,x = 1 — Yhx^Y -Px,w 

Note that if r = 0, i.e. there are no light particles, then the two-species ASEP is 

simply the usual ASEP. 

3.3. The Matrix Ansatz. We now explain the Matrix Ansatz, a technique intro¬ 
duced in [DEHP93] for computing the stationary distribution of the ASEP. We also 
present Uchiyama’s generalization of the Matrix Ansatz l.'cliOS to the two-species 
case. Finally we will recall the solution to the Matrix Ansatz which was found in 
|USW04j . 

For convenience, we now set u = 1. Derrida, Evans, Hakim, and Pasquier 
[ DEHP93 ] proved the following theorem. 

Theorem 3.3. [ DEHP93] Suppose that there are matrices D and E, and vectors 
{W\ and \V) such that the following relations hold: 

• (W\(aE-~/D) = (W | 

. (PD-5E)\V) = \V) 

• DE — qED = D + E 

Let Zx = (W\(D + E) N \V). Then in the ASEP on a lattice of N sites, the steady 
state probability of state (t\, , Tx) is equal to 

wn;h[»(n = *)p+i(n = °i£]|y') 

Z N 

For example, the steady state probability of state (•, o, o, •, •) is equal to ( W \ DE ^ DD \ V ) 
Uchiyama generalized the Matrix Ansatz to the setting of the two-species exclusion 
process. 

Theorem 3.4. [Uch08] Suppose that there are matrices D, E, A, and vectors (IT| 
and \V) such that the following relations hold: 

• (WKaE-'yD) = (W \ 

. U3D-6E)\V) = \V) 

• DE — qED = D + E 

• DA = qAD + A 

• AE = qEA + A. 
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Let ZN tr = [y r ] ^ W ^ D (^A^\v) Then in the two-species ASEP on a lattice of N 
sites with precisely r light particles, the steady state probability of state (ti, ..., r n ) 
is equal to 

{W\ nilitlfo = 2 )D + Ifo = 1)A + 1 (Ti = 0)E]\V) 

EN,r 

For example, the steady state probability of state (2, 0,1,1, 2, 2) is equal to \deaadd\x) 


3.4. A solution to the Matrix Ansatz. Now consider the following tridiagonal 
matrices, which were introduced by Uchiyama, Sasamoto and Wadati in [USW04| FI 


d 


" d ^ 
a 0 

4 

1 

o 


r e 11 

e o 

4 

0 • 


> 
u 0 

A 

A 


e o 

4 

e l 


• o 

d\ 

4 

, e = 

0 

e i 

4 ‘ 











where 


4 := <fl n (a,b,c,d) 


4 : = e n( a , b, c, d) 


(1 — q 2n ~ 2 abcd)( 1 — q 2n abcd ) 

x \bdfa + c) + (b + d)q — abcdfb + d)q n ~ 1 — {bd(a + c) + abcdfb + d)}q 
— bd(a + c)q n+1 + ab 2 cd 2 (a + c)q 2n ~ 1 + abcdfb + d)q 2n ], 
q n ~ 1 

(1 — q 2n ~ 2 abcd){ 1 — q 2n abcd ) 

x [ac (6 + d) + (a + c)g — abcd(a + c)g n_1 — {ac(b + d) + abcd(a + c)}q 
— ac(b + d)g n+1 + a 2 bc 2 dfb + d)q 2n ~ 1 + abcd{a + c)g 2n ], 


4 : = 4( a > 6 > c > d ) = !> 


4 := e*^(a, 6, c, d) = -g n ac, 


4 := 4(°Ac,d) 
4 := 4(°Ac,d) 


(fbd 

[1 — q n ac)[l — q n bd) 

(l - q n ac)(l - q n bd) An ' ^ 


A -72 . A n (n, 6, c, d) 

(1 — g n_1 a&cd)(l — q n+l )( 1 — q ,n afe)(l — g n ac)(l — q n ad) (1 — q n bc)( 1 — g n fed)(l — q n cd ) 
(1 — g 2n_1 a6cd)(l — q 2n abcd) 2 ( 1 — q 2n+1 abcd ) 


1 Actually we have slightly modified the definitions of djj, dj, ejj, and e■ which appeared in [ USW04; 

but it is easy to check that both the matrices we are using here and the original matrices of |USW04| 
satisfy Lemma 13.51 
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Lemma 13.51 below uses the following change of variables between a, b, c, d and 


( 4 ) 

( 5 ) 

( 6 ) 

( 7 ) 


a 

c 

b 

d 


l-g-a + 7 + a/(1 - g - a + 7) 2 + 4ay 
2 a 

1 — q — a + 7 — a /(1 — q — a + 7) 2 + 4«7 
2 a 

1 - q - p + 5 + y/(l - q - p + S) 2 + 4/35 
2 P 

1 - q - (3 + 5 - y/{l - q - /3 + 5 ) 2 + 4/35 

2/3 


Note that this change of variables can be inverted via 


( 8 ) 

(9) 


1-9 

a = --, 

1 + etc + a + c 

— (1 — q)ac 

7 = tz - r — 

1 -\- clc -\- d H - c 


1 H - bd H - b d 
— (1 — 

1 -J - bd b ~jr d 


Lemma 3.5. [USW04] Let D = y4-(l + d) and E = + e )? where 1 is the 

identity matrix. Also, use the equations through (J7J) to express D and E in terms 
of a , /3, 7 , S, and q. Let (W| = (1, 0, 0,...) and \V) = {W\ T . Then D, E, (W\, and 
\V) give a solution to the Matrix Ansatz of Theorem, 1,9.,91 


Moreover, Uchiyama |Uch08] observed that if A := DE — ED , and D and E satisfy 
DE — qED — D + E, then it follows that DA = qAD + A and AE = qEA + A. 
Therefore we have the following. 


Corollary 3.6. LetD,E, (W | , and\V) be as in Lemma \3.<% and set A DE—ED. 


Then D, E, A, (W |, and \ V) give a solution to the Matrix Ansatz of Theorem\3.f 


In the rest of the paper, we will use the solution of the matrix ansatz in Corol¬ 
lary 13.61 which in turn fixes our definition of the following versions of the partition 
function. (Note that one gets a definition of the partition function from every solu¬ 
tion to the Matrix Ansatz; using a different solution will simply rescale the partition 
function.) 


Definition 3.7. Let Z N (ff) = (W\(£D+E) N \V) be the fugacity partition function for 
the ASEP where D, E, {W\, and \V) are as in Lemma liTR We also let = 

[y r ] ^ ^ ^ e the fugacity partition function for the two-species ASEP. 

Finally we set Z N = Z N ( 1) and Z N>r = Z Nr ( 1) and refer to these quantities as 
partition functions. 


4. Koornwinder moments and partial Motzkin paths 

We start by reviewing the connection between the partition function Z N of the 
ASEP and moments of Askey-Wilson polynomials. We will then explain how to 
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interpret Zn as a generating function for weighted Motzkin paths. Finally we will 
use the celebrated lemma of Karlin-McGregor-Lindstrom-Gessel-Viennot [ KM59b . 
IKM59allLm78llGV85 on determinants and non-intersecting lattice paths to express 
Koornwinder moments in terms of partial weighted Motzkin paths. 


4.1. The partition function of the ASEP, Askey-Wilson moments, and 
Koornwinder moments. Using results from [ IJSW04 J. we proved in [ GSSW12] 
that there is a close relationship between the fugacity partition function of the 
ASEP and Askey-Wilson moments /i at. 


Theorem 4.1. |CSSW12[ Theorem 1.11] Recall that Z N (f] a, (3,^/, 8; q) = (W\(fD + 
E) N \V) be the fugacity partition function of the ASEP, where D, E, fW\, and \V) 
are as in Lemma \3 .51 Then the N th Askey-Wilson moment is equal to 


(1 - q) N 

p N (a,b,c,d\q) = Z N y 


2 N i N 


-l;o,/3,7,5;g),0 


where i 2 = — 1 and 

( 10 ) 

i -q 


a = 


1 — ac + ai + ci 


7,0 = 


i - q 


1 — bd — bi — di 


7, 7 = 


(1 — q)ac 


1 — ac + ai + ci 


7,8 = 


(1 — q)bd 


1 — bd — bi — di 


Recall from Lemma 12.71 our determinantal expression for Koornwinder moments. 
Because of the close relationship between the partition function of the ASEP and 
Askey-Wilson moments (see Theorem 14. II) . we define another version of Koornwinder 
moments as follows. 


Definition 4.2. Given a partition X = (Ai, A 2 ,..., X m ), we define the Koornwinder 
moment at q = t to be 


( 11 ) 


KxiO 


det(Z\ 

i+m—i+m—j (0)%=i 

det(Z 2m _ i _ i (0)™=i 


where Z N is as in Definition \ d. 7\ 


The fact that we refer to the quantities AR(£) as Koornwinder moments is justified 
by the following result. 


Proposition 4.3. Let A = (Ai,..., A m ) be a partition. We have that 


M\(a, b, c, d\q) 



|A| 

Ji A (-l;a,/3,7,5;g), 


where |A| = Yhi an d a, (3, 7 , 8 are related to a, b, c, d as in fflOl) . 


2 The formula we give here does not include the factor {afi — , because we are working 

here with the partition function coining from the D and E from Lemrna l3.51 rather than the rescaled 
version which was used in ICSSW12l . 
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Proof. Let 9 = So Zjv(— 1; a, (3, 7, 8]q) — 9 N hn((i, b, c, d|g). Then we have that 

det(Z Ai+m _ i+m _ J -(-l))^ =1 


2 i 

Kx( 


- 1 ) = 


det(Z 2m -i-j( 1 ))ij=i 
det (Q-Pi+m-i+m-j) 


H'Xi+m—i+m—j )i,j=l 


det( 0 -( 2 m --i ) /i2m _ i _ j )- =1 

_ g- (Al+ - + A " I + m(m - 1) ) det (^ +m -i + m-j)Z=l 
e - mim - l ) det ( /i2m _ i _ i )™ =1 

= g|—(A iH -hA^) m 


□ 


From now on, when we refer to Koornwinder moments, we will be referring to 
Definition 14.21 

Because of its probabilistic interpretation, the quantity Zn must be positive when¬ 
ever we specialize the parameters a, (3,^,8, and q to be positive numbers between 0 
and 1. Moreover, we proved in [jCWll :, [CW13] that (up to a normalizing factor) Zn 
is a polynomial in ct, f3, 7 , 5, q with positive coefficients; it can be expressed as a sum 
over staircase tableaux. This prompted Rains to conjecture the following. 

Conjecture 4.4. The Koornwinder moment K\(ff) is a polynomial in a, (3, 7 , 8, q, f 
with positive coefficients (up to a normalizing factor). 

Since Theorem 11.11 implies that Rqjv-r,o,o,...,o)(£) is proportional to the fugacity 
partition function Z^ r {f() for the two-species ASEP, it follows that the Koornwinder 
moment -Rpv-r, 0 , 0 ,...,o) is positive when we specialize the parameters a, (3, 7 , 8, q to be 
positive numbers between 0 and 1. Moreover, in jCMWj . we will give a tableaux 
interpretation of Z/v, r , which implies that Zjv,r(£) is a polynomial in £,a,/ 3 , 7 , 8, 
and q, with positive coefficients. Finally, see Section [8] for some evidence towards 
Conjecture 14.41 when q — £ — 1. 

4.2. Weighted Motzkin paths and the partition function of the ASEP. 

Recall that a Motzkin path of length A is a path in the xy plane from (0, 0) to (A, 0) 
which consists of steps northeast (1,1), east (1,0), and southeast (1, —1), and never 
dips below the :r-axis. One often associates a weight to each of the three kinds of 
steps, based on the height at which the step begins. The weight of a given Motzkin 
path is then the product of the weights of all of its steps, and such a path is called a 
weighted Motzkin path. If C = {c l3 ) is a tridiagonal matrix rows and columns indexed 
by the non-negative integers, then we can use each nonzero entry c i3 to weight a step 
in a Motzkin path that starts at height i and ends at height j. We call such a Motzkin 
path a C-Motzkin path , though we will sometimes drop the C if it is understood. 

The following lemma is obvious. 

Lemma 4.5. Let C be a tridiagonal matrix as above, and let (W\ = (1,0,0,...) 
and \V) = (W\ T . Let Z N = ( W\C n \V). Then Z N is the generating function for all 
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C-Motzkin paths of length N, i.e. it is the sum of the weights of all C-Motkin paths 
from (0, 0) to (. N , 0). 


We will now give a general result relating a certain ratio of determinants to partial 
Motzkin paths, which we will use to subsequently relate Koornwinder moments to 
partial Motzkin paths. Note that it is well-known that there is a link between (par¬ 
tial) Motzkin paths and moments of (one-variable) orthogonal polynomials, see for 
example jViel Chapter 1, Proposition 17] and also |Vie85j . However, our treatment 
here will be self-contained. 

We say a partial Motzkin path is a path in the xy plane from (0, 0) to (N, r) which 
consists of steps northeast (1,1), east (1,0), and southeast (1,-1), and never dips 
below the x-axis. As before, we can use entries cy,- of a tridiagonal matrix C to 
associate the weight cy,- to a step which starts at height i and ends at height j. The 
weight of a partial Motzkin path is then the product of the weights of all of its steps. 
We have the following result. 

We define the following ratio of determinants: 


( 12 ) 



Z\ 1 +2n-2 

-2’Ai+2n—3 • 

-2^1+71— 1 

det 

'2'A2+2n—3 

-2’A2+2n—4 • 

Z\ 2 + n —2 



£\n+n— 2 • 

Z \n 



Z 2n—2 

Z 2n—3 ■ 

• Z n -1 


det 

Z 2n—3 

Z 2n—A ■ 

Z n-2 



Z n -1 

Z n -2 

■ Z o 



Theorem 4.6. Let C be a tridiagonal matrix, and let C Motz(lV, r) be the generating 
function for all C- partial Motzkin paths which start at (0,0) and end at ( N,r ). Let 
(W | = (1, 0, 0,...), \V) = ( W\ T , and \V r ) = (0,..., 0,1, 0, 0,... ) T , where the 1 is 
in the rth position (and coordinates are indexed by the non-negative integers). Let 
Z N = {W\C N \V) and k r = UiZo c M+1 . Then we have that 


(13) 


fc Motz(lV, r) = fl}V\C N \V r ) = K (JV _ r , 0 ,o.o), 

r\j r 


where there are precisely r 0 ’s in (N — r, 0, 0,..., 0). 
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Proof. The leftmost equality in (jT3j) is obvious. To relate these quantities to the 
ratio of determinants 


det 


(14) 


IC(N- 


(N-r, 0,0,...,0) 


-2 N+r 

Z-N+i —1 • • 

Z N 

Z 2r -1 

Z 2r -2 

■ Z r -1 

Z r 

Z r -1 

•• tq 


det 


z 2r 

Z 2r —i .. 

Z r 

txi 

to 

* 

1 

Z 2r - 2 ■ ■ 

■ Z r _i 

•• tq 

^ .. 

1 

• ^0 


at the right, we now use the well-known Karlin-McGregor-Lindstrom-Gessel-Viennot 
Lemma (KM59bl [KM59al. !Lin731 1GV85 J (which we will henceforth refer to as the 
KMLGV Lemma). This lemma will give a combinatorial interpretation of both the 
numerator and denominator. 


(r,r) 



(r,r) 



Figure 2. An acyclic directed graph (we consider all edges to be 
directed either east, northeast, or southeast) corresponding to the de¬ 
nominator of (U3l) when r = 4, and the unique collection of pairwise 
vertex-disjoint paths from {A 0 ,..., A r } to {B 0 ,..., B r } in this graph. 


Consider the acyclic directed graph shown at the left of Figure [2] (for the case 
that r = 4). Although we haven’t shown the orientations of the edges, we consider 
all edges to be directed either east, northeast, or southeast. The points A 0 , A lr .., 
A r = B r , F> r _ i, L? r _ 2, ..., Bq are at the lattice points (0, 0), (1,0), ..., (r, 0),(r + 
1, 0),(r+ 2, 0), ..., (2r, 0), respectively. We weight each edge of the graph which goes 
from height i to height j by the entry Cy of C. Then we define an (r + 1) x (r + 1) 
weight matrix M den = (M dcn ) as follows. Its rows and columns are indexed by 
{0,1 ,... , r}, and M den is defined to be the weights of all paths in the graph from 
Ai to Bj. In other words, M den is the generating function for all C-Motzkin paths 
of length 2 r — i — j. Applying Lemma [4.51 we have that M den = and hence 

M den is the matrix appearing in the denominator of 

But now by the KMLGV Lemma, the determinant det M den is the generating func¬ 
tion for pairwise vertex-disjoint path collections from {A 0 ,..., A r } to {B 0 ,..., B r }. 
It is easy to see that there is only one such path collection, which is shown at the 
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right of Figure [2] Note that [Kra991 Theorem 11] can also be used to give a formula 
for this Hankel determinant. 



(r,r) 



Figure 3. An acyclic directed graph (we consider all edges to be 
directed either east, northeast, or southeast) corresponding to the nu¬ 
merator of (TT3T) when r = 4. Now there are multiple collections of 
pairwise vertex-disjoint paths from {A 0 ,..., A r } to {B 0 ,..., B r }, but 
all of them use the set of paths from Ai to B L for 1 < i < 4 shown at 
the right of the figure. 


Now consider the acyclic directed graph shown at the left of Figure [3] (for the case 
that r = 4). Again we consider all edges to be directed either east, northeast, or 
southeast. The points A 0 , A lr .., A r = B rj B r _ i, B r _ 2 , ..., B 0 are at the lattice 
points (—N + r, 0), (1, 0), ..., (r, 0),(r + 1, 0),(r + 2, 0), ..., (2r, 0), respectively. We 
again weight edges of the graph using entries of C. Then we define the (r + 1) x (r + 1) 
matrix M num = (M™ m ), with rows and columns indexed by {0,1 ,..., r}, as follows. 
Mff im is defined to be the weights of all Motzkin paths from Ai to Bj. Using Lemma 
14.51 it follows that M num is the matrix appearing in the numerator of (TT3l) . 

By the KMLGV Lemma, the determinant det M num is the generating function 
for pairwise vertex-disjoint path collections from {A 0 ,..., A r } to { B 0 ,..., B r }. It 
is easy to see that such path collections must consist of the paths from Ai to B{ 
for 1 < i < r shown at the right of Figure |3l together with an arbitrary Motzkin 
path from Aq to Bq which goes through the point (r, r ) and whose last r steps are 
southeast. 

We now have combinatorial interpretations for both det M num and det M den . Com¬ 
paring them, we see that the ratio is equal to the generating function for 

the partial Motzkin paths from (—N + r, 0) to (r, r) divided by the weight of the 
unique partial Motzkin path from (0, 0) to (r, r). Clearly the generating function for 
the partial Motzkin paths from (—N + r, 0) to (r, r) equals the generating function 
for the partial Motzkin paths from (0,0) to ( N,r ). And the weights of the up steps 
in our Motzkin paths are c^j+i, so the weight of that unique partial Motzkin path is 
Co, 101,2 • • • c r —i^ r . This completes the proof. □ 

Corollary 4.7. Now we choose the tridiagonal matrix C = £,D + E, where D and E 
are as in Lemma \S. ,51 Let CMotz(A^, r) be the generating function for all C- partial 
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Motzkin paths which start at (0, 0) and end at (. N , r ). Then we have that 
K{ N - r ,o,o,...,o) = ^(W\C N \V r ) = jpC Motz(iV, r), 

r\jy* !\/rp 

where there are precisely r 0 ’s in (N — r, 0, 0,..., 0), and k r = n[=o (£ — 

To prove Theorem 11.11 we now need to relate (W\(fD + E) N \V r ) to the fngacity 
partition function for the 2-species ASEP. 


5. A Jacobi-Trudi formula for general Koornwinder moments 


In this section we will prove a general Jacobi-Trudi type result in Theorem 15.11 
which expresses the quantity IC\ (see (1121) ) in terms of the quantities /C( m ,o,...,o)- Our 
proof techniques build on those used in the proof of Theorem 14.71 We again let 
C = (cij) be a tridiagonal matrix. Edges of graphs that we will construct here are all 
directed northeast, east, or southeast, and each edge from height i to j is weighted 
by dj. 

Theorem 5.1. Let A = (Ai,..., A n ) be a partition. Then 

(15) JC\ = det(/C(A i +j-j ) o,o,...,o))i(j=n 
where the term on the right-hand-side has precisely n — j 0 ’s. 

Note that Theorem 15.11 immediately implies Corollary 15.21 which expresses a gen¬ 
eral Koornwinder moment K\(C.) in terms of the “complete homogeneous” Koorn¬ 
winder moments AT( m)0 ,. ,.,o)( 0 - 

Corollary 5.2. Let A = (Ai,..., A n ) be a partition. Then 

(16) im) = det(A' (A , +3 _ j ,„.„ l 0 ) (e))?j_ 1 , 

where the term on the right-hand-side has precisely n — j 0 ’s. 


We now prove Theorem 15.11 

Proof. We start by analyzing the left-hand side of (fT5l) . By definition, we have that 


det 


(17) 


K 


(Al, —, A n) 


-^Ai —|— —2 

A^Ai~t-2 n —3 • • 

Z^Ai+n—1 

^A2+2n—3 

£a 2+ 2„-4 

^\2+n— 2 

£\n+n—l 

-ZAn+n—2 



det 


-Z^2n—2 

^271-3 • • 

• %n- 1 

^2n—3 

^271—4 • • 

Z n -2 

Zn-1 

2 

. Z o 


Using the same arguments as in the proof of Theorem 14.71 we interpret the matrix 
in the numerator of (TT7T) as the weight matrix associated to the acyclic directed graph 
at the left of Figure 01 where edges from height i to height j are weighted by c ir Note 
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*g+l 

<- 1 - 

*2+2 

*4+3 



Figure 4. An acyclic directed graph (all edges are directed either 
east, northeast, or southeast), whose weight matrix is given by (TTTlh 
when n — 4. 


that the points A 1; ..., A n , and B n , ..., Bi all lie on the x-axis, with the B^s at 
unit distance apart. The positions of the Aj’s have been chosen so that the distance 
between Aj and B n is \ -\-n — i. By Lemma [4.51 the ij entry Z\. +2n -i-j of the matrix 
in the numerator of (1T7l) is equal to the generating function for Motzkin paths from 
Ai to Bj , hence this matrix is the weight matrix associated to the graph at the left 
of Figure |U By the KMLGV Lemma, its determinant is the generating function for 
pairwise vertex-disjoint path collections from {Ai, ..., A r } to {£>i, ■ • •, B r }. Note 
that all such path collections must use the southeast edges which are shown in bold 
at the left of Figure |4j Moreover, as we showed in the proof of Theorem 14.71 the 
determinant in the denominator of (TT7|) is equal to the weight of the path collection 
shown at the right of Figure [4] Therefore the ratio of the determinants is equal to 
the generating function for pairwise-disjoint non-intersecting paths from A 1: ..., A n 
to the lattice points Bi,..., B 4 in Figure [5j divided by the product of the weights of 
all the up steps in the path collection shown at the right of Figure [4j That product 
is equal to k\k 2 ■ ■ ■ k n _ i, where k r = ni^c M+ i. 



*3+1 

<- 

*2+2 

*4+3 


Figure 5. An acyclic directed graph, shown for n = 4. Edges are 
oriented northeast, east, or southeast. 
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We now consider the right-hand side of (fi~5lh We have that 
(18) 



^(Ai,0,...,0) 

A+Ai+1,0,...,0) 

^d(Ai+n—1) 

det(/C(A i +j-i,o,o,...,o))”_/=i — det 

^C(a 2 -i,o,...,o) 

/C(A2,o,...,o) 

^d(A2+n—2) 


fc(\n-n+ 1,0,...,0) ^C(An-n+2,0,...,0) • 

£(*») 


where the partitions in column j of the matrix contain precisely n — j 0’s. Using 
Theorem 14.71 we have that (1181) is equal to the determinant of 

-A— CMotz(Xi + n — 1, n — 1) -A— CMotz(X\ + n — 1, n — 2) ... CMotziXi + n — 1, 0) 
-A— CMotz(X 2 + n — 2 , n — 1 ) -A— CMotz(X 2 + n — 2 , n — 2 ) ... CMotz(X 2 + n — 2 , 0 ) 

En — 1 l^n —2 

-A—CMotz{X n , n — 1 ) -jA—CMotz{X n , n — 2 ) ... CMotz( A n , 0 ) 

which is equal to 
(19) 


CMotz( Ai +n — l, n — 1) CMotz(Xi + n — 1, n — 2 ) 
CMotz( X 2 + n — 2 , n — 1) CMotz( X 2 + n — 2,n — 2) 


CMotz(Xi + 7 i 
CMotz( X 2 + 7 i 


1 , 0 ) 

2 , 0 ) 


n—1 

n 


^=1 


i 

ke' 


CMotz(X n ,n — 1) CMotz(X n ,n — 2 ) 


CMotz(X n , 0) 


But now it’s clear that the matrix in (1191) is the weight matrix for the directed graph 
shown in Figure EH so we are done. □ 


6. From partial Motzkin paths to the two-species ASEP 


Our goal in this section is to prove the following. 


Theorem 6.1. Define D, E, (W\, and \V) as in Lemma lTR and set A = DE — ED. 
Using the change of variables from (jSJ) and (f9j) , define 


( 20 ) p r = 

Then we have that 


(i -g) r 

IKte-tfac) 


(21) (W\(fD + E) N \V r )-p r 


(1 -q) r a r (l-g) r 
k r ni: 0 1 «+? , 7)' 

{W\(£D + E + yA) N \V) 
[V 1 {W\A r \V) 


In other words, (W\(fD + E) N \V r ) ■ p r is the partition function for the two-species 
ASEP on a lattice of N sites with precisely r light particles. 


Note that our main theorem, Theorem 11.11 is an immediate consequence of Theo¬ 
rem 14.71 and Theorem 16.11 

We will actually prove a refinement of Theorem 16.11 To state the refinement, we 
need to introduce some notation. Given a word X in {D, E } N , let S r (X) be the set 
of words which can be obtained from A" by replacing precisely r letters in X by an 
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A. Given a word Z e S r (X), let D(Z) (respectively, E(Z)) be the set of positions 
of letters that were D (respectively, E) in X and became A in Z. Let 

in v E (Z)= ^ \{i£D(Z)UE(Z)\i<j}\. 


j€E{Z) 


Theorem 6.2. Let p r = a r {l — q) r . For any word X in {D, E} N , we have that 

(Z)„,\D(Z)UE(Z)\ ( W\Z\V) 

(W\A r \V)' 


( 22 ) 




'r 


(W\X\V r )-p r = ^2 
zes r {x) 

We next explain how Theorem 16.21 refines Theorem 16.11 
Lemma 6.3. Theorem 1 6. 2\ implies Theorem \6.1[ 


Proof. Define \X\ D to be the number of D' s that occur in a word X. We take f|22]) 
and sum it over all words X in {D, E} N , keeping track of the number of ZTs in X, 
and obtaining 
(23) 


E (} x '°(W\X\V T )-~p r 

X&{D,E} N 



E f 


\X\ Dq iuv E (Z) a \D(Z)\^\E(Z)\ 


xe{D,E} N zeSr(x) 


(W\Z\V) 

(w\A r \vy 


On the left-hand side of f[23|) . we obtain (W\(£D + E) N \V r ) ■ p r = (W\(fD + 

W-ft-lCK + W- 

Now we will analyze the right-hand side of ([23]) . Let us fix a word Z e {D, E , A\ N 
which contains precisely r A’s, and determine the coefficient of ^y\A^v) • Note that 
\X\ D = \Z\ D + D{Z). Also, will appear 2 r times on the right-hand side - 

this is the number of ways that each of those A’s in Z could have come from a D 
or an E in a word X G {D,E} N . Now note that the statistics E(Z), \D(Z)\, \E(Z)\ 
only depend on the substring of r A’s in Z , and which letter each of those A’s 
originally corresponded to in X. Let us consider each of the A’s in Z from left 
to right. If the first A came from a D (respectively, an E), then we pick up a 
factor of af (respectively, 7 ) in f\ x \oq mVE ( z ')a\ D( ' Z ' > \ r y\ EI ' Z) \. If the second A came 
from a D (respectively, an E), then we pick up a factor of af (respectively, gy) 
in qinv E (z) a \D(z)\^\E(z)\ _ ^.j ie q comes from the inversion between the first and 
second A. And in general, if the ith A came from a D (respectively, an E), then we 
pick up a factor of af (respectively, g'^y). Therefore the coefficient of 011 

the right-hand side of (1231) is (af + y)(a£ + gy) • • • (ct£ + g r ^ 1 y). It follows that the 
right-hand side equals [y r ] • 11 [=0 + 3 * 7 )- 

Comparing our expressions for the left-hand side and right-hand side yields (j2Tj) . 

□ 


Theorem 6.4. To prove Theorem \6.A it suffices to prove Theorem. \6.S\ in the case 
that X = D N . In other words, it suffices to prove that 


(24) 


(W\D N \V T ) . Pr 


v (W\Z\V) 
2- (W\A*\V) 

Z&Sr(D N ) ' ' 


(W\(D + yaA) N \V) 
[V 1 {W\A r \V) 
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Note that the second equality of (1241) is obvious. Theorem 16.41 is a consequence of 
the following two lemmas. 

Lemma 6.5. If Theorem \6.2\ is true for all words X of the form E £ D m , then it is 
true for all words X in the letters D and E. 


Proof. Let U be an arbitrary word in D and E. By Corollary 13.61 the matrices I). 
E, and A satisfy the relations of Theorem 13.41 In particular we have that DE = 
qED + D + E. Note that by repeatedly applying the relation DE = qED + D + E 
(to replace instances of DE by qED + D + E), we can put U into “normal form” - 
that is, we can write U as a linear combination of words of the form E e D m . We will 
prove Lemma 16.51 by induction on the number of times one must apply the relation 
DE = qED + D + E to put U into normal form. The base case is the words which 
already have the form E e D m . 

For the inductive step, let us write U = XDEY , where X and Y are words in 
D and E. We have that XDEY = qXEDY + XDY + X'EY. By the induction 
hypothesis, Theorem 16.21 is true for each of the words XEDY, XDY, and XEY. 
We want to show that (l22li also holds for XDEY. 

Let us first analyze the right-hand side of (l22l) when applied to the word U = 
XDEY. We will write each word Z G S r (XDEY ) as ZxDEZy or ZxAEZy or 
Z x DAZy or ZxAAZy , where Zx and Z Y have been obtained from X and Y, re¬ 
spectively, by replacing some of the letters by A’s. We will also write ZxZy to 
denote the word obtained from one of the four words above by deleting the two 
letters between Z x and Z Y . 

We have that Y^zes r (u) g mVg a I 1 7 1 I ^ is equal to 


E 


Z=Z x DEZ Y eS r (U) 


inv E (Z x Z Y ) \D{Z X Z Y )\ \E(Z x Z y )\ (W\Z\V) 

q J (W\A r \V) 


+ E 

Z=Z x AEZ Y £S r (U) 

+ E 9 

Z=Z x DAZ Y £S r (U) 

+ E 9' 

Z=Z x AAZ Y eSr{U) 


n \E(Z Y )\+imr E (Z x Z Y ) \D(Z x Z y )\+1 \E(Z x Z Y )\ mz\v) 

1 1 (W\A*\V) 


\E{Z Y )\+\D(Z x )l+\E(Z x )\+inv B (Z x Z Y ) \D(Z X Z Y )\ \E{Z X Z Y )\+1 (W\Z\V) 

1 {W\A r \V) 


2\E(Z Y )\+\D(Z x )\+\E(Z x )\+l+inv F AZxZ Y ) \P(Z X Z Y )\+1 ,\E(Z X Z Y )\+1 ( W \Z\V) 

1 (W\A r \V)' 


Note that in the first three terms above, where Z has the form ZxDEZy, ZxAEZy, 
and ZxDAZy, we can apply the relations DE = qED + D + E, AE = qEA + A, and 
DA = qAD + A, to rewrite (W\Z\V) as (W\Z x (qED + D+E)Z Y \V), (W\Z x (qEA+ 
A)Zy \V), and (W\Zx(qAD + A)Z Y \V), respectively. 

Now let us analyze the left-hand side of (122|) when applied to the word U = 
XDEY. Applying the re-writing rule DE = qED + D + E and using the induction 
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hypothesis, we have that (W\X(DE)Y\V r ) ■ p r is equal to 


(q(W\X(ED)Y\V r ) + (W\XDY\V r ) + (W\XEY\V r )) ■ p r 


E * 

ZeSr(XEDY) 

+ E 

ZGSr(XEY) 


inv E (Z)+l \D(Z)\ \E{Z)\ ( W\Z\V) 

1 (W\A r \V) 

inv E (Z) \D(Z)\ \E(Z)\ (W\Z\V) 

q 1 {W\A r \V) 


+ E 

ZGSr(XDY) 


Yv E (Z) \D(Z)\ \E(Z)\ ( W \ Z \ V ) 

7 {W\A r \V) 


E 


Z=Z x EDZ Y eS r (XEDY) 


n inv E (Z x Z Y )+l \D(Z X Z Y )\ \E(Z X Z Y )\ mz\v) 

‘ ’ (W\A'\V) 


+ 


E 


Z=Z x EAZ Y eS r (XEDY) 


n \E(Z Y )\+l+inv E (Z x Z Y ) \D(Z x Z y )\+1 \E(Z x Z y )\ (W\Z\V) 

q 7 (W\A r \V) 


\D(Z x )\+\E(Zx)\+\E(Z Y )\+l+inv B {ZxZ Y )JD(Z x Z Y )\ n/ \E(Z x Z Y )\+l mm 

(1 W\A r \V) 


7 


+ E 9 

Z=Z x ADZ Y eS r (XEDY ) 

I V |D(Z x )|+|E(Z x )|+2|E(Zy)|+l+inv B (Z x Zy) |D(Z x Z y )|+l |B(Z x Zy)|+l 

^ ^ 7 (Wl/W) 

Z=Z x AAZ Y eS r (XEDY) \ \ \ / 


+ 


E 


Z=Z x DZ Y £S r {XDY) 


+ 


E 


Z=Z x AZ Y &S r (XDY) 


inv E (Z x Z Y ) \D(Z X Z Y )\ \E(Z x Z y )\ (gj 
|B(Z x )|+m v B (Z x Z x ) |D(Z x Zy)|+l |E(Z x Zy)| 


+ E 

Z=Z X EZ Y eS r (XEY) 

+ E 9 

Z=Z x AZ Y eSr(XEY) 


uiv E (Z x Z Y )\D(Z x Z Y )\ \E(Z x Z y )\ mz\v) 

1 7 {W\A r \V) 


\D(Z X )\+\E (Z x )\+\E(Z Y )\+inv E (Z X Z Y ) \D(Z X Z Y )\ \E (Z x Z Y )\+1 jg) 

' (W' r | J 4 r |y) 


Now if we compare our expressions for the left-hand side and right-hand side of 
(122]) . we see that they are equal. This completes the proof. □ 

Lemma 6 . 6 . If Theorem \6.2\ is true for all words X of the form D N , then it is true 
for all words X of the form E f D m . 


Proof. We will prove that Theorem 16.21 is true for any word U containing exactly I 
E’s by induction on I and on n, the number of applications of DE = qED + D + E 
which are necessary to put the word in normal form. 

If U has the form X(DE)Y, then the proof of Lemma 16.51 shows that the theorem 
holds for U, by rewriting U = X(DE)Y = qX(ED)Y + XDY + XEY. (Note that 
the process of rewriting does not increase the number of E ’s in the word.) 
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Otherwise we can assume that U has the form U = E e D m . Let us first analyze 
the left-hand side of (1221) for the word U. By Corollary 13.61 and Theorem 13.41 we 
have the relation (W\E = ±(W\ + %{W\D. 

Using the induction hypothesis, we have that (W\E e D m \V r ) ■ p r is equal to 


{-(W\E e ~ 1 D m \V r ) + ^(W\DE e - 1 D m \V r )) ■ ~p r 
a a 


E 


mv E (Z) a \D(Z)\-l \E(Z)\ 


Z&SriEl-'D™) 


C W\z\v) 

(W\A r \V) 


+ 

+ 


E 

Z=DZ x £S r (DE e - 1 D m ) 


E 

Z=AZ x £S r (DE l - 1 D m ) 


Q 

q 


inv E (Z x ) \D(Z X )\-1 \E(Z X )\+1 (W\Z\V) 

1 (W\A r \V) 

\E{Z x )\+inv B (Z) \D(Z X )\ \E(Z X )\+1 mz\v) 

' (W\Ar\V)' 


Analyzing the right-hand side of (I22|) for the word U , and again using the relation 
(W\E = ±{W | + 1 (W\D, we have that 


E 


Z=EZ x eS r (E e D m ) 


mv E (Z x ) \D(Z X )\ \E(Z X )\ ( W \ Z \ V ) 

q 7 (W\A r \V) 


+ 


E 


Z=AZ x eSr(E t D™) 


\E(Z x )\+mv E (Z x ) \D(Z X )\ \E(Z X )\+1 jwqzjD 

1 ’ (W\A'\V) 


E 

Z=EZ x &S r (E t D m ) 


V x )„Wz x) u E{ z x) \ l a m Z x\V) + imDZxW) 

1 (W\A r \V) 


q mVE ^'a 


+ 


E 


Z=AZ x &S r (E l D m ) 


\E(Z x )\+inv E (Z x ) \D(Z X )\ \E{Z X )\+1 mz\v) 

1 ’ (w\A'\vy 


Comparing our expressions for the left-hand side and right-hand side of 
see that they are equal. This completes the proof. 


we 

□ 


Theorem 6.7. The equation (124|) holds. More specifically, 


(W\D N \V r )-p r 


( W\(D + yaA) N \V) 
[V 1 (W\A r \V) 


Section [T] will be devoted to the proof of Theorem 16.71 Note that once we have 
proved Theorem 16.71 this theorem together with Theorem 16.41 will imply Theorem 
16.21 And then Theorem 16.21 and Lemma [6.31 will imply that Theorem 16.11 holds. 


7. The proof of Theorem 16.71 

In this section we will prove Theorem 16. 71 We will primarily work with the matrices 
d and e defined in Section 13.41 so we begin by writing down the relations satisfied by 
these matrices. 
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Lemma 7.1. Let d and e be defined as in Section \3J\ let D, E, {W |, and \V) be 
defined as in Lemma HOI let A = DE — ED; and let \V r ) = (0,..., 0,1, 0, 0,... ) T , 
where the 1 is in the rth position, so that |V°) = \V). Then we have the following 
relations. 


d A = qAd 
Ae = qeA 

de = ged + (1 — q) 1 
d k e = q k ed k + (1 - g^d^" 1 
d\V) = (b + d)\V) - bde\V) 

(W\e= (a + c)(W\-ac(W\d 

d|U) = 4_ 1 \v t ~ 1 ) + 4\v“) + 4\v i+1 ). 


Proof. This is a straightforward exercise; most relations follow immediately from the 
corresponding relations that D and E satisfy. □ 


Let [k] ■= 1 + q + • • • + q k 1 be the g-analogue of the number k, and let 
[NW- 

[r] q \[N-r] q \- 


N 

r 


9 


Theorem 7.2. We have that 


(25) 


(W\d N \V r ) 


(W\(d + yA) N \V) 
[y J (W\A r \V) 


N 1 (W\A r d N ~ r \V) 
r \ g (W\A r \V) 


Lemma 7.3. Theorem\7.2\ implies Theorem\6 . 7 
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Proof. Suppose that Theorem 17.21 is true. Then 


(W\D N \V r ) • p r = (W\ ^ + d l N N \V r ) • a r { 1 - qy 


0--q) 


= a 


N / \ 


i =0 
N 


-n^ N\,rAW\(i + yA)'\V) 


= a r{i -«r A X , m 


i =0 
r-N N 


{W\A r \V) 


o r (i -q) 

(W\A*\V) j^\i 


E (^)blW[(1 -q)D- 1 + VA]'\V) 
% 0 1(1 ~ q)D ~ 1 + vAm 

a r ( 1 - q) r ~ N 


(W\A r \V) 


■m<wm~q)D+yAf\v) 


a 


= W ] 


{W\A r \V) 

(.D + yaA) N \V) 


[y r ](W\(D + yA) N \V) 


(W\A r \V) 


□ 


To prove Theorem 17.21 we will find an explicit formula for the quantity in (j25j) (see 

Corollarv l7.10p . and use recurrences to show that both (hh|d iV |f/ r ) and 

are equal to this explicit formula. Note that the second equality in Theorem 17.21 fol- 
lows by repeated application of the relation d A = qAd. 


N 

r 


(W / |A r d JV-r | V) 
(W\A r \V) 


7.1. A formula for ( W\A r d N r \V). First we will prove a recurrence for (W\A r d m \V) 
(where we are thinking of m = N — r). 

Proposition 7.4. We have that 

(w\A r i'"\v) = 6 + 1 : mA^- 2 w). 
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Proof. In what follows, we will use the relations from Lemma 17.11 Note that 

(W\A r <T\V) = (W\A r 6 m ~ l ((b + d)\V) - bde\V)) 

= {b + d){W\A r d m - 1 \V) - bd{W\A r d m ~ l e\V) 

= {b + d)(W\A r d m ~ 1 \V) - bd (q m - 1 (W\A r ed m - 1 \V) - ( q m ~ l - l)(W\AA m ~ 2 \V)) 

= {b + d)(W\A r d m ~ 1 \V) - bdq m ~ 1 (W\A r ed m ~ 1 \V) + bd{q m ~ 1 - l){W\A r 6 m ~ 2 \V) 

= {bp d){W\A r d m ~ 1 \V) - bdq m+r ~ 1 (W\eA r 6 m ~ l \V) + bd{q m - 1 - l)(W\A r d m - 2 \V) 

= {b + d)(W\A r d m - 1 \V) - bdq m+r - 1 [(a + c){W\A r d m ~ l \V) - ac{W\dA r d m - l \V)] 

+ bd{q m ~ 1 - l){W\A r d m ~ 2 \V) 

= {b + d){W\A r d m - 1 \V) - bd(a + c)g m+r - 1 (PL|A r d m - 1 |I/) + abcdq m+2r - l {W\A r d rn \V) 
+ bd{q m ~ 1 - 1){W\A V d m ~ 2 \V). 


□ 

Our next goal is to get an explicit formula for (see Corollary I7.10p . 

Towards this end, we make the following definition. 

Definition 7.5. Define F m (y ) by F 0 (y ) = l, F m (y) = 0 if m < 0, and 

(26) F m {y) = {b + d-y{a + C )q m - 1 )F m . 1 {y)) + {q m - 1 -l){bd-acq m - 2 y 2 )F m _ 2 {y) 

for m > 0. 

Corollary 7.6. We have that 

{W\A r d m \V) _ F m (bdq r ) 

(W\A'\V) ~ n?=~ 0 \l-abcdq^y 

Proof. Corollary 17.61 follows from Proposition 17.41 and Definition 17.51 by comparing 
the recurrences and base cases. □ 


m 

i 


Theorem 7.7. Set 

m 

B m {b,d) = J2 

i =0 L 

Then for m > 0, we have 

m 

(27) F m (y) = Y, 1 .- 1 ) 1 


b l d m 1 and A m (a,c ) = ^ 


»=o 


m 

i 


a l c m ~\ 


V? 


i =0 


m 

i 


q^y l A i (a,c)B rn _ i (b,d). 


Before proving Theorem 17.71 we first state a few useful lemmas. Lemma 17.81 is a 
simple exercise (and is well-known). 
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Lemma 7.8. 

(28) 

(29) 

(30) (1 - O 
Lemma 7.9. We have that 


m 
i 

m 
i 

m — 1 

i 


= q‘ 


m — 1 

i 

m — 1 

i 


= (1 ~q m ~ l ) 


+ q' n 

+ 

9 

m 

i 


m — 1 
i — 1 


m — 1 
i — 1 


- 1 9 


(b + d)B m (b , d) = d) + (1 — q m )bdB m _i{b , d), and 

(a + c)A m (a, c) = An+i(a, c) + (1 - g _m )arn4 m _i(a, c ). 

Proof. It suffices to prove the first equation. We will use Lemma [7.81 Note that 


(b + d)B m (b,d) = J2 


m 

i 


m +1 


E 


6 i d m+1_i 

i=0 L J 9 i=l L 

m 


m 

i — 1 


i jm+l-i 


b l d 


- 1 9 


d m+ i + 6 m+1 + 

Z=1 

m 

d m+ 1 + 6 m+1 + 

m+l 

E 


Z=1 


+ 


m 
i 

- 1 9 

m + l 
i 


m 

i — 1 


b i d m+1 - i 


- 1 9> 


+ (1 - g m “* +1 ) 


m 

i — 1 


L 

9/ 


^d m+1 - 


Z— 0 L 

m+1 

E 

z=0 
m+l r 

E 

i=0 


m + l 
i 


6 i d m+1_i 






m + l 
i 

m + l 

i 


j =i 

m—1 


m 


j -1 


Vd m+1 - j 


ijm+l—i 


b'd 


+ E (1-9”+ 


z=0 


m—1 r 


6*d m+1_ * + (1 - q m ) 


i =o 


m 


m — 1 


- 1 9 

6* +1 d m - i 


6 i+1 d m - i 


= B m+1 (6, d) + (1 - g m )6d5 m _i(&, d). 


□ 


We now prove Theorem 17.71 

Proof. The case m = 0 is trivial. We prove the proposition by induction. Suppose 
that the proposition is true for all n < m. Then we have the following (see the end 
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of the list of equations for justifications of each line): 


('b + d - y(a + c)q m )F m _ 1 (y) 


m— 1 


2—0 


m — 1 

i 


(■b + d-y(a + c)q m *) J^(-l) 

m—1 

£<-i> 

2—0 
m—1 

2=0 
m 

E<-!) 


m — 1 

i 


q ( 2 ) y l A (a, c ) (6, d) 


- 1 9 


i=0 


l 

m—1 


i 

m — 1 


+E<-i) 

2 

m 

- Ewe "*- 1 

i=0 

-E(- 1 )V n_I 

m — 1 


m — 1 

i 


i 

m — 2 


E(-d 

i=0 

+E*- 1 ) 

2 

m 

i=i 

2 

m 

E(-d 

o 

E<-i) 

2 


q^)y'Ai(a, c)(b + d)B Tn -i-i{b 1 </) 

"* ~ 1 | 9 (y +m - 1 9 i+1 (a + c)A(a, c)B m _,_,(6, d) 
q < 3y l A i (a,c)B m _ i (b,d) 
q^y'Ma, c)(l - g”*- i - 1 )6dS m _ i _ 2 (6, d) 
g( 2 )r/ l+1 ^l m (a, c)5 m _ i _ 1 (6, d) 
m . 1 g(a)|/* +1 ac(l - g“*)A-i(a, d) 

X 

J 9 

q^)y l A i (a,c)B m _ i (b,d) 
q&tfMa, c)(l - d) 

g( J 2 )y J A j (a,c)B m _ j (b,d ) 

m q(^y l+1 ac(l - q rn ~ 1 )A i _ 1 (a,c)B, m _ i _ 1 (b,d) 


m — 1 

:i -1 


1=0 


m — 1 

i 

m — 2 

i 


i — 1 

+ Q r 


m—1 
i — 1 


q^y l Ai(a,c)B, m _i(b,d) 


- 1 


g(*)j/%(a,c)(l - g”*- 1 )6dS m _ i _ 2 (6,d) 

<? 

m — 2 
i — 1 


g(2V +1 ac(l - g m d) 
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i =0 


m 


q^)y l A i (a,c)B m _ i (b,d) 


m — 2 

i 


+(i-nwi(-i) 

i 


q^ 2 )y l Ai(a, c)B m ^ 2 {b, d) 
g 

m — 2 


J 


A(°> 2 (b, d) 


Z(-1)‘ 


i =0 


m 


g(»)j/%(a, c)S m _,(6, d) + (1 - q m ~ x ){bd - acg— 1 2 y 2 )F m _ 2 (y) 


- 1 <? 


For the first equality, we used the induction hypothesis, and for the second equality 
we simply distributed terms. For the third equality we used Lemma 17.91 For the 
fourth equality, we changed the index of summation in the third line, and we also 
used Lemma 17.81 


q^-i-i) 

m — 1 

= (i — q m ~ l ) 

m — 2 


i 

9 

i 


and 


(1 - Q l ) 


m — 1 

i 


= (i - q 


m—l\ 


m — 2 
i - 1 


- 1 9 


For the fifth equality, we combined two terms. For the sixth equality, we used 
Lemma 17.91 and we changed the index of summation in the third line. For the 
seventh equality, we used the induction hypothesis. 

Now we have shown that 


(b + d - y(a + c)q m 1 )F m _ 1 (y) 


is equal to 


EM)' 

1=0 


m 

i 


q($y l Ai(a, c)B m _i(b , d) + (1« q m x ){bd - acq m 2 y 2 )F m _ 2 (y ). 


Comparing this to the defining recurrence for F m (y) (Definition I7.5|) . we have proved 
the desired explicit formula for F m (y). □ 


Corollary 7.10. We have that 


{W\A r d m \V) 


Er= 0 (-i ) 1 


m 


g(a) (bdq r yAi(a, c)B m _i(b, d) 


(W\Ar\V) n™o \l~abcdq 2r + i ) 

Proof. This follows from Corollary 17.61 and Theorem 17.71 


□ 
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7.2. A formula for (W\d N \V r ). First we will prove a recurrence for (VF|d Ar |’F r ). 
Towards this end, we state a simple lemma, which follows from the definitions of the 
tridiagonal matrices d and e. 

Lemma 7.11. 


d\V r ) = (1 - q 2r ~ 1 abcd)\V r ~ 1 ) - bdq r e\V r ) + R r \V r ), 

where 

T — 1 

R r = - 2 r -2 (?(& + d) — q r ^ (b + d)abcd + bd(a + c)(l — q r )) . 

Theorem 7.12. The quantity fW\d N \V r ) satisfies the following recurrence. 

(1 - abcdq r+N ~ 1 )(W\d N \V r ) = (1 - q 2r ~ 1 abcd)(W\d N ~ 1 \V r ~ 1 ) - bdq r ( 1 - g JV ” 1 )(lT|d Ar - 2 |^ r ) 

+ {R r - bdq r+N ~\a + c))(lT|d iV - 1 |l/ r ). 

Proof. In what follows, we will use Lemma 17.111 and Lemma 17.11 - in particular the 
relations d fc e = q k ed k + (1 — q k )d k ~ l and (W|e = (a + c)(W\ — ac{W\d. 


(W\d N \V r ) =(1 - q 2r - 1 abcd)(W\d N - 1 \V r ~ 1 ) - bdq r (W\d N ~ 1 e\V r ) + J R r (IT|d JV - 1 |1/ r ') 

=(1 - g 2r_1 a6cd)(lT|d JV_1 |ld r_1 ) + R r {W\d N ~ 1 \V r ) 

- bdq r+N ~ 1 (W\ed N ~ 1 \V r ) - bdq r ( 1 - g Ar - 1 )(IT|d iV - 2 |1/ r ) 

= (1 - q 2r - 1 abcd)(W\d N - 1 \V r ~ 1 ) + R r (W\d N ~ 1 \V r ) - bdq r+N ~\a + c)(W\d N - 1 \V r ) 
+ abcdq r+N ~ 1 (W\d N \V r ) - bdq r ( 1 - q N - 1 )(W\d N ~ 2 \V r ). 

Collecting like terms now yields the recurrence in Theorem 17.121 □ 


Theorem 7.13. We have that 
(W\d m+r \V r ) = 


m + r 
r 


Fm{bdq r 


n^o^l -abcdq 2r+i )' 


Note that once we have proved Theorem 17.131 combining this theorem with Corol¬ 
lary 17.61 will prove Theorem 17.21 and hence Theorem 16.71 


Proposition 7.14. Define 

C(m,r) = --. 

nS (l - “bcd q ^‘) 

Then C(m,r ) satisfies the recurrence 

(1 —abcdq 2r+m ~ 1 )C(7n, r ) = Afim, r)C(m, r—l)—A 2 (m, r)C(m— 2 , r)+A 3 (m, r)C(m— 1 , r), 
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where 

(1 — q 2r ~ l abcd)(l — q r ) 

1 ~~ (1 - q m + r ) 

bdq r {l — g m-1 )(l — q m ) 

2 ~ (1 - q m+r ) 

(— bd(a + c)q m ~ 1+2r + i? r )(l — q m ) 

— -;-;— - . 

(1 - q m+r ) 

Equivalently, F m (bdq r ) satisfies the following recurrence. 

(1 - q m+r )( 1 - abcdq 2r ~ 2 ) F m (bdq r ) = (1 - <f )(1 - abcdq 2r+rtl - 2 )F m (bdq r ~ 3 ) 

- bdq r { 1 - g m “ 1 )( 1 - q m )(l - abcdq 2r ~ 2 ){l - a6crfg 2r+m - 2 )F m _ 2 (Mg r ) 
+ (1 - q m )(-bd(a + c)q 2r+m - 1 + R r )( 1 - abcdq 2r - 2 )F m _ 1 {bdq r ). 


Lemma 7.15. 


To prove Theorem 7.13, it suffices to prove Proposition 7. If 


Proof. To prove Theorem 17.131 we need to show that 


(31) 


(w\a m+T \v T ) 


m + r 
r 


C(m, r). 
q 


It is easy to verify this equation for m + r < 1. Moreover it is straightforward to 
verify that the recurrences in Proposition 17.141 are equivalent to the recurrence for 
(W\d m+r \V r ) (see Theorem 17.121) . provided that we have (I3B . □ 


Therefore to complete the proof of Theorem 17.131 we need to prove Proposition 
17.141 In particular, we will prove the recurrence for F m (bdq r ), which follows from 
Proposition 17.161 below when y = q r . 

Proposition 7.16. F m (bdy ) satisfies the following recurrence. 

(1 - yq m )( 1 - y 2 q~ 2 abed)F m (bdy) = (1 - y 2 q m ~ 2 abcd)( 1 - y)F m (bdy/q) 

— bdy( 1 — q m )( 1 — g m_1 )( 1 — abedy 2 q m ~ 2 )( 1 — abcdy 2 q~ 2 )F m _ 2 (bdy) 

+ (1 — q m ){y [b + d + q^ 1 bd(a + c)] — y 2 [g^ 1 (l + q m )bd(a + c) + q~ 2 (b + d)abcd ] 
+ y A \q m ~ 3 ab 2 cd 2 (a + c)] }F m -i(bdy). 

We will prove Proposition 17.161 by taking the coefficient of y n in the recurrence 
and showing that it is an identity. Towards this end, we define 

(32) X min = (-1) 

which is the coefficient of y n in F m {bdy). 


q(^b n d n A n (a, c)B m _ n (b, d), 
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Lemma 7.17. Taking the coefficient of y n in the recurrence of Proposition \ 7.16 
yields the following equation. 

(33) 

X m ,n q 0>bcdX m>n — 2 U X mn _i -f- q abcdAm^n -3 

=q~ n X m , n - q~ n+1 AVi - q m ~ n abcdX m , n - 2 + q rn - n+l a.bcdX m . n _ 3 
-(1 - q m )( 1 - q m ~ 1 )bdX m - 2 ,n- 1 + q~\ 1 - g m )( 1 - q^ab 2 cd 2 X m _ 2 , n _ 3 
+ (1 - q m ){ 1 - q m ~ l ) \q m ~ 2 ab 2 cd 2 X m _ 2jn _ 3 - q m ~ 4 (abcd) 2 bdX m _ 2>n _ 5 ] 

+ (1 — q m ){b + d)X m -i >n -i + q *(1 — q" l )bd{a + c)X m _ 

-(1 - g m ) [g _1 (l + q' n )bd(a + c) + q~ 2 abcd{b + <f)] A m _i )n _ 2 

+ (1 — q m )q m 5 ac(a + c)( 6 d) 2 A m _i ) 7 l _ 4 . 

Proof. The proof is straightforward so we omit it. □ 

Note that the total degree in a and c in X m ^ n is n, while the total degree in b and 
d in X mjH is m + n. Moreover, each term in (1551) has one of the following properties: 

• the total degree in a and c is n, and the total degree in b and d is m + n, or 

• the total degree in a and cis n — 1, and the total degree in b and d is m+n — 1. 
So we can split (1551) into two equations based on which of these two conditions holds. 
We get equations (155)) and (155]) . respectively. 

(34) q 2 {l - q~ n )X m ^ n - abcd( 1 - q m ~ n+2 )X mtn _ 2 - q( 1 - q m )bd(a + c)X m _ hn _ 1 

+ (1 — q m )abcd{b + d) X m _i n _ 2 = 0. 

(35) 

(q- n+1 ~ q m )X m , n -1 + (q m - 2 ~ q m ~ n+1 )abcdX m , n _ 3 - (1 - q m )(b + d)X m _ ljU _ 1 
+ g -1 (l - q m )( 1 + q m )bd{a + c)X m _ 1<n _ 2 - g m_3 ( 1 - q m )ac(a + c)(bd) 2 X m _ hn _ 4 
+ (1 - q m )( 1 - g m - 1 )6d{X m _ 2 , n _ 1 - (q- 2 + q m - 2 )abcdX m _ 2 , n _ 3 

+ q m -\abcd) 2 X m _ 2 , n _ 5 } = 0. 

Now it is clear that Proposition 17.161 follows from Lemmas 17.181 and 17.191 below. 
To prove the lemmas, we will take the coefficient of a^b n+l c n ~^d m ~ l in each of (154)) 
and (1551) . obtaining identities involving g-binomial coefficients which can be proved 
by hand or with any computer algebra system. 

Let x(m,n,i,j) denote the coefficient of a^b n+l c n ~^d m+n ~ l in X mjn (q). It is easy 
to check that 

(36) x(m, n, i, j) — (—l) n g( 2 ) 


m 


m — n 


n 

n 


i 


3 


Lemma 7.18. Equation (154)) is an identity. 
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Proof. It suffices to prove that taking the coefficient of a-ffi" + *c n J d m 1 in (134|) yields 
an identity. When we take this coefficient we obtain the following equation. 

(1 — q~ n )x(m,7i, i,j) — (g -2 — q m ~ n )x{7n, n — 2, i + 1 , j — 1 ) 

— g _1 (l — q m ) [x(m — 1 ,n — 1, i,j ) + x(m — 1 , n — 1 , i, j — 1)] 

+ g” 2 (l — q m ) [x(m — 1 , n — 2, i + 1 , j — 1) + x(m — 1, n — 2, i, j — 1)] = 0. 

Using (1361) now yields the equation 


1/9 




777 


m — 77 


n 

n 


i 


_ j 

L 

9 


9 



-2 


777 


777 — np 2 


77 — 2 

n — 2 

9 

i + 1 

9 

_ j - 1 . 


J !/<? 


-1 


+ g U 2 -j- 2 (l-g m ) 


(1 - g m ) 


m 


n 


m 


71 


m — 7i 

i 

m <— n + 1 

i 



71—1 

+ 

77 — 1 


u 

. j - 1 . 

1/9 

3 

1 / 9 / 


+ 


777, — 71 + 1 

i + 1 


77—2 

3~ 1 


= 0. 


1/9 


Simplifying this, we get the equation 

^- 1(1 _ g m-n+2 )(1 _ q3)(l _ q n- 3) ^ 


q 2n -\l-q~ n ) 


+ q 2n -\l - q~ j )(l - q~ n+j ) 


(1 - g i+1 )(l - qm-n+l-i^ 

1 1 


+ 


q- *(2 -q- j -q- n+j ) 
= 0, 


]_ _ q m—n+l—i ^ _ gl+* 

which is trivially true. □ 

Lemma 7.19. Equation (j35|) is an identity. 

Proof. We use the same strategy as in the proof of Lemma 17.181 Now we take the 
coefficient of b n ~ 1+l c 71-1 ^ d m ~ z in (j35lh obtaining the following equation. 

(q~ n+l — q m )x(m, n — 1, i,j) + (g m ~ 2 — q m ~ n+1 )x{7n, n — 3, i + 1, j — 1) 

+ (1 — g m )(l — q m ~ 1 ){x(m — 2 , n — 1 , i — 1 , j) — (q~ 2 + q m ~ 2 )x(m — 2, n — 3, i,j — 1) 

+ q m ~ A x(m — 2 , n — 5, i + 1, j — 2)} — (1 — q m ){x(m — 1,77— 1 , i,j ) 


+ x(m — 1 ,77 — 1, i — 1 ,j) 


1 + q r 

Q 


[x(m — 1,77 — 2, i,j — 1) + x(m — 1,77 — 2, i,j )] 


+ q m 3 [x[m — 1, n — 4, i + 1, j — 1) + x(m — 1, n — 4, i + 1, j — 2)]} = 0 
Again it is straightforward to show that this is an identity. □ 

8. KOORNWINDER MOMENTS AT £ = q — 1. 

Our goal in this section is to show that when one evaluates Koornwinder moments 
at £ = q = 1, one obtains a beautiful multiplicative formula in terms of hook lengths. 
This formula is given in Theorem 18.51 below. 
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To make sense of this specialization, however, we need to be careful. It is clear from 
Section [XT that the Uchiyama-Sasamoto-Wadati solution to the Matrix Ansatz has 
poles when q — 1, which is related to the fact that some of our formulas for Koorn- 
winder moments - in particular Theorem 11.11 - are only valid for q ^ 1. However, 
there are other solutions to the Matrix Ansatz which are well-defined at q = 1. We 
will review one such solution here, use it to define the partition function and Koorn- 
winder moments (which are the same as our previous definitions, up to a global scalar 
factor), and then state and prove Theorem 18.51 

In this section we will abbreviate K\(l\ a, /3, 7 , 5; 1) by writing K\. 


8.1. The Matrix Ansatz at q — 1. We give now a solution to the Matrix Ansatz 
at q — 1, which comes from [ USW041 Section 5]. We define tridiagonal matrices T> 

El El 0 

El E } E\ 

0 E\ E\ ’ 


|y) — ( 1 , 0 , 0 , • • • ) T , 

ph _ Z 3 + 7 + + 2^7 + 7<S) 

(a+ 7)00 + 5) 

E i = — J— l( n + l)(^ + n )} 1/2 , 

£ '» = ^[(« + i)U + «)] 1/2 , 

with 

_ a + /S + 7 + h 
(a + 7)(/3 + 5) ' 

We now use the above solution to the Matrix Ansatz to define the partition function 
and Koornwinder moments at q = £ = 1. 

We define the following notation: 

a = Q + 7)(/ 3 + 4) 

a(3 — 7 <f 

Definition 8.1. When q = 1, we define the partition function of the ASEP by 
Z N = S N (W\(V + ^) Ar |V)@ And we define the Koornwinder moment K\ by fTTT]l . 
after setting £ = 1 and using the above definition of Z^. 

3 We have included the constant S N in the definition of the partition function so that this 
definition of Zn equals the specialization of the fugacity partition function from Definition 

13.71 at £ = q = 1. 


and E, and vectors (W\ and \V) by 

£>5 D l 0 •••' 

Dl D\ D[ 


(37) 


V = 


where 


ZT 

n 

D i 

n 

D b = 


0 D\ D\ 


(W\ = ( 1 , 0 , 0 , • • •), 


a + 5 + n(af3 + 2 a5 + 7 5) 


a 


a + 7 

5 

JTd 


(o + 7)(/ 3 + h) 

[fn + l)(x + n)] 1/2 , 

[(n + l)(x + n)] l/2 , 
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Let C = ( Cij ) be the tridiagonal matrix defined by 
1 1 

CU = S(D\ + E\) = S(x + 2n) 
cy-i = S 2 (A-, + = S 2 i(x - 1 + i). 

Lemma 8.2. We have that Zn = fW\C N \V). 

Proof. Note that for every step from height i to height % + 1 in a Motzkin path, there 
must be a corresponding step from height * + 1 to height i in the Motzkin path. It 
follows that (W\C N \V) = S N (W\{V + S) N \V). □ 

Proposition 8.3. We have that 

/1\t\ N ~f i /N\ N .~. 1 

K (N-r, 0,0 . 0 ) = sN ~ r y r J = ( aj g - 7 6)N-r (r ) II («+/ 5 +7+<5+*(«+7)(^+<J))> 

where there are precisely r 0 ’s in (N — r, 0, 0,..., 0). 

Proof. If we apply Theorem 14.61 and use the fact that each c^+i = 1, we get 
(38) C Motz(iV, r ) = (W\C N \V r ) = K (N - r ,o,o,...,o). 

But now note that there is a simple recurrence for partial Motzkin paths: 

C Motz(iV, r) = C Motz(A r — 1, r)c r>r + C Motz(iV — 1, r — l)c r _i jr + C Motz(iV — 1, r + l)c r+ i 
= CMotz(iV — l,r)S(x + 2 r) + CMotz(iV — l,r — 1) + 

CMotz(iV — 1 ,r + 1 )S 2 (r + l)(x + r). 

This implies the result. □ 

Lemma 8.4. We have that ifyv-r,o,...,o) = Zn, t , where the partition (N — r, 0,..., 0) 
has precisely r 0 ’s. 

Proof. When q — 1 we have that A = DE — ED = D + E. Therefore 

[; y r }(W\{yA + D + E) N \V) [y r ]{W\(y + 1) N (D + E) N \V) { N r )Z N 

N ’ r {W\A r \V) (W\A r \V) Z r ' 

Since Z N = K N (by (ITT]) ). Proposition l8.3l tells us that Z N = ^ a+ ^J+^N a+ ' ) ' >i ' l3+S ' > ■ 

Therefore Z N r = ( r ) W=r . Combining this with Proposition 18.31 

we see that K(N-r,o,...,o) = Zn, t - D 

8.2. Koornwinder moments at q — £ — 1. We identify each partition A = 

(Ai,...,A n ) with its Young diagram, namely a left-justified array of cells with n 
rows and A, cells in the ith row for 1 < % < n. For each cell z of the Young diagram 
in position the hook H\(z) is the set of cells (a, b) such that a — i and b > j, 

or a > i and b = j. The hook length h\{z) is the cardinality of H\(z). 

Theorem 18.51 is the main result of this section. 
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Theorem 8.5. For any partition A = (Ai,..., A n ), we have 

n_1 ” (x + Ai - Aj + j - i - 1)(A* - Aj + j - i) 


k x =s | a ' n<-+ Kz) - 1 ) n n 

z£ A i =1 j=i +1 


(x + j 


1 )U 


To prove Theorem 18.51 we need a few lemmas. 

Lemma 8.6. Let A = (Ai,..., A m ) and v = (Ai,..., A m , 0). Then 

(x + Aj + m — i) (Aj + m — i + 1) 

2=1 


k v =K\ . n 


(x + i — l)(i) 


Proof. By Corollary 15.21 and Lemma 18.41 we have that 

K v = det 

where 

( ,\ Vi+m-i 

T ) n 

m + 1 — 7 J 

J ' i=m+l-j 


Note that 


( \ . | -I \ A i+m—i 

x TXr) n <*+<> 

J ' o =m +l-j 


for i rn + 1, and 


M m+1 , = S j ~ m+1 


l m+l,j 

Therefore 

Meanwhile we have that 
where 


0 A J 1 if j — m + 1 

m + 1 — j) otherwise. 

\m+1 


det(M«)’”t 1 1 = det(M«)”> =1 . 
/U = det (MD3 = i, 

• \ Ai+m—2—1 


77‘) n (*+ 

J 7 £=m—j 


Clearly MA = ' M ir Therefore 

(m + 1 — j) (x + m — j) 


K\ = det 


(Aj + m + 1 — i ) (x + A, + m — i) 


■Mi 


u 


i,j =1 


nr=i(™+ 1 - i)7+ m - j) 


nr=i(A. + m + 1 - i)(x + Aj + m — i) 
n;=i(77+j-i) 


rii = i(7 + m + 1 - i)(x + Aj + m — i) 
This proves the lemma. 


-detfMy)™., 

K„. 
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Lemma 8.7. Let X = (Ai,..., A n ) and is — (Ai + 1, ..., A n + 1). Then 

n 

K v = K\S n • + Xi + 


n — i . 


i=i 


Proof. By Corollary 15.21 and Lemma 18.41 we have that 

K u = det (M i3 )l j=v 

where 

( \ i i i *\ Ai+n— i 

„ , ) n <*+*>■ 

s i=n—j 


Meanwhile we have that 


K x = det (Mr)l j=v 


where 


/ \ \ Xi+n-i -1 

ML = + ” - A n (* + Q- 

U d J e=n _j 


We can now write 


n-j 7 V n- j - 1 


\i~\~n—i 


n 


Z=n-j 


= + A, + n - i) + MM d + A ‘ + n ~; . 

3 l ’ 3+ x + n-j-1 

But now since the determinant is alternating in the columns of the matrix 


we have that 


K v = det(My) = det (SMhfx + A i + n — i ))™ j= i 

n 


i= 1 


S n K x + Aj + n - i). 


i=1 


This proves the lemma. 


□ 


Proof of Theorem \8.5[ We start with a partition A = (Ai,..., A„). If n — 0 then 
K\ = 1 and we are done. Otherwise if X n = 0, we apply Lemma [8761 and get 

n—1 


^A -^(Ai,...,A n —l) ’ | J[ 


= K 


(Al ,...,An—l) 


1=1 

n—1 

n 

i =1 


(x + Aj + n - 1 — i ) (A* + n — i ) 

(x + i- 1 )(i) 

(x + Xi — X n + n — 1 — i) (A i — X n + n — i) 
(x + n — i — l)(n — i) 
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By induction, we know that 


n —2 7i—l 


k { n 

zeA *=i j=*+i 

and therefore 

ri—1 n 

Kx= s |A| n<*+ Hz) -1) n n 

zGA i= 1 j=i+l 

Finally if X n > 0, we apply Lemma 18.71 


(x + \-Xj+j-i - 1)(A i - A j + j - i) 


(■x + j-i - 1 ){j -i ) 

(x + \i — \j + j — i — l)(Aj - A j + j — i) 


(x + j 


1 )U 


K\ — Ai-i,...,An-i) ' S n ' + Aj + n — i). 


2=1 


It is easy to see that 


I^[(x +/r(^) - 1) = n( x + n “ 0 • n (x + h(z)~ 1). 

2GA j=1 zg(Ai— l,...,A n —i) 

By induction, we know that the theorem holds for ^(Ai-i^.^An-r), an d therefore we 
get the desired result for K\. □ 

We now explain why Conjecture 14.41 is true when £ = q — 1. First note that by 
iterating Lemmas 18.61 and 18.71 we can prove the following result. 

Proposition 8.8. Let v = (Ai+1,..., A m +1, 0,..., 0) and A = (Ai,..., A m , 0,..., 0), 
where there are r 0’s at the end of each partition. Then 

^ ^ T-r Aj + m + r + 1 — i . , 

K u — K\ ■ —-— , _ , 1 - —(x + A i + m + r -i). 


2=1 


A* + m + 1 


Corollary 8.9. Conjecture 4-4 is t rue when £ = q = 1, namely the Koornwinder 
moments at this specialization are polynomials in a, (3,^,5 with positive coefficients. 


Proof. Note that we can repeatedly apply Proposition 18.81 to an arbitrary partition 
(Ai,..., A„), so as to express K\ in terms of K(o,o,...,o)- It follows from the definition 
that A(o,o,...,o) is equal to 1 , and applying the proposition preserves the property of 
being a polynomial with positive coefficients. 

□ 
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